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Summary 


Static  strength  and  damage  tolerance  analysis  methods  for  skin  modifications  to 
antenna  installations  on  the  fuselage  of  commuter  airplanes  are  described  in  this 
report.  The  methods  are  implemented  in  the  analysis  program  of  RAPIDC 
Version  1 .0  for  commuters.  Three  types  of  antenna  installation  in  terms  of  the 
shapes  of  mounting  plate  are  first  defined.  The  applicability  of  antenna  systems 
within  the  analysis  capability  of  RAPIDC  is  next  addressed.  Static  strength  and 
damage  tolerance  analysis  procedures  are  then  described.  The  damage 
tolerance  analysis  procedure  is  demonstrated  through  an  illustrative  example. 

1.0  RAPIDC  Capabilities 

The  analysis  of  a  modified  fuselage  skin  to  antenna  installation  in  RAPIDC  for 
commuter  airplanes  is  an  addition  to  the  RAPID  program  originally  developed  for 
common  repairs  to  large  transport  aircraft.  The  analysis  capability  of  RAPIDC 
Version  1.0  for  antenna  installations  is  limited  to  the  types  and  systems 
described  below. 


1.1  Types  of  Antenna  Installation 


Three  common  types  of  antenna  installation  are  considered.  They  differ  only  in 
the  shapes  of  mounting  plates:  rectangular,  circular,  and  elliptical  shape  as 
shown  in  Figure  1 .  The  mounting  plates  are  mechanically  fastened  to  the  skin 
around  the  cutout. 


(a)  Rectangular 


(b)  Circular 


(c)  Elliptical 


Figure  1 .  Three  types  of  antenna  installation 

In  RAPIDC  Version  1.0,  only  the  circular  and  rectangular  installations  are 
considered.  The  elliptical  installations  will  be  included  in  the  next  version 
release. 


1.2  Antenna  System  Applicability 

RAPIDC  Version  1.-0  excludes  the  case  of  skin  bending  that  is  induced  by 
antenna  vibrations  subjected  to  aerodynamic  loading.  Therefore,  communication 
antennas  (typically  5”  long  by  3”  wide  by  10”  tall)  are  excluded.  However, 
antenna  systems  such  as  ADF  (12”  long  by  6”  wide  by  1”  tall),  GPS  (4”  long  by  4” 


wide  by  0.5”  tall),  Transponder  and  DME  (5”  long  by  2”  wide  by  4"  tall),  and. 
TCAS  are  covered  in  the  analysis. 

2.0  Static  Strength  Analysis 

Typical  antenna  installations  on  the  fuselage  of  commuter  airplanes  generally 
involve  cutting  a  hole  in  the  skin  for  the  antenna  connector.  A  mounting  plate  is 
then  mechanically  fastened  to  the  skin  around  the  skin  cutout.  To  assess  the 
static  strength  of  the  modified  skin,  three  independent  criteria  are  used  to 
evaluate  the  margins  of  safety  of  the  mounting  plate  and  fasteners:  the  mounting 
plate  allowable,  the  joint  allowable,  and  the  shear  allowable. 

2.1  The  Mounting  Plate  Allowable  and  Margin  of  Safety 

To  assess  the  loss  of  skin  strength  due  to  antenna  connector  hole  cutout,  the 
margin  of  safety  of  the  mounting  plate  is  calculated  as 

Pu 

Margin  of  Safety  =  ~  - 1 

where  Ppu,  the  mounting  plate  load  allowable,  is  calculated  using  the  equation 


Ppu-  F,u  (Wp  D)  tp 

in  which  Ftu,  wp,  D,  and  tp  are  the  ultimate  tensile  strength,  width,  cutout 
diameter,  and  thickness  of  the  mounting  plate,  respectively. 

The  skin  internal  force,  Pp  ,  is  calculated  using  the  equation 


Pp  =  au  D  ts 

where  <ru  and  ts  are  the  design  ultimate  tensile  stress  and  thickness  of  the  skin, 
respectively. 

2.2  The  Fastener  Joint  Allowable  and  Margin  of  Safety 

In  the  vicinity  of  an  antenna  installation  on  the  skin,  skin  stresses  are  transferred 
to  the  mounting  plate  through  fasteners  via  fastener  shear  and  hole  bearing. 
The  fastener  together  with  the  skin  and  mounting  plate  represents  a  fastener 
joint.  A  joint  can  only  transfer  a  certain  amount  of  load  until  it  fails.  The  point  at 
which  this  fastener  joint  fails  is  the  joint  allowable.  The  calculation  of  the  joint 
allowable  is  based  upon  two  loading  conditions,  fastener  shear  and  hole  bearing. 
Two  hole  types  also  have  to  be  considered  when  determining  the  joint  allowable, 
the  straight  shank  and  countersunk. 
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2.2.1  Fastener  Joint  Allowable  in  Straight  Shank  Holes 

Straight  shank  holes  are  used  for  protruding  head  fasteners  and  for  the  mounting 
plate  that  have  a  flush  head  fastener  installed  but  are  not  countersunk.  The  joint 
allowable  is  the  lower  of  the  shear  allowable  or  the  hole  bearing  allowable.  An 
allowable  is  calculated  for  the  skin  and  mounting  plate  that  the  fastener  goes 
through. 

The  single  shear  allowable  for  straight  shank  holes  is  calculated  using  the 
following  equation 


PSu  -f"su  x  Af  XSCF1 

The  mounting  plate  shear  allowable  for  straight  shank  holes  is  calculated  using 
the  following  equation 


Psu=2x Fsu  X  Af  X S CF2 

In  both  equations,  Fsu  is  the  ultimate  shear  strength  of  the  fastener  material,  Af  is 
the  cross-sectional  area  of  each  fastener,  Scfi  is  the  single  shear  correction 
factor  and  Scf2  is  the  double  shear  correction  factor.  Scfi  and  Scf2  are  used 
only  for  solid  rivets  and  can  be  found  in  MIL-Handbook  5F  Table  8.1 .2.1  (b). 

The  hole  bearing  allowable  for  the  straight  shank  hole  is  calculated  using  the 
following  equation 


Pbru=Fbruxdxt 

in  which  F^  is  the  ultimate  bearing  stress  of  the  plate  (skin  or  mounting  plate) 
material,  d  is  the  fastener  hole  diameter,  and  t  is  the  thickness.  Currently,  the 
Fbru  is  for  the  case  E/D  (edge  distance  to  hole  diameter)  equal  to  2.0. 

The  joint  allowable  for  a  given  joint  is  the  lower  of  Psu  or  Pt,rudenoted  by  Pjoint. 

2.2.2  Fastener  Joint  Allowable  in  Countersunk  Holes 

The  fastener  joint  allowables  for  countersunk  holes  are  different  from  the  straight 
shank  holes.  They  are  determined  by  tests  and  can  be  found  in  MIL-Handbook 
5F  Section  8. 

2.2.3  Fastener  Joint  Margin  of  Safety 

The  fastener  joint  margin  of  safety  is  determined  by  summing  each  fastener  joint 
allowable  in  the  skin  and  mounting  plate,  and  this  is  done  for  each  side  of  the 
mounting  plate. 


3 


^total  2  \  j joint  )n 
n»1 


in  which  Ptotai  is  the  total  fastener  joint  load,  (Pjoint)n  is  the  fastener  joint  load  for 
the  nth  fastener  and  K  is  the  number  of  fasteners. 

The  Ptotai  for  the  skin  and  mounting  are  then  compared.  The  smaller  of  the  two  is 
the  fastener  joint  allowable  for  that  side  of  the  mounting  plate  and  is  used  in 
determining  the  margin  of  safety  for  the  fastener  joints. 

An  applied  load  is  needed  to  determine  a  margin  of  safety.  That  load  is  the 
ultimate  applied  load  to  the  structure,  or  if  that  is  not  known,  then  the  tensile 
ultimate  strength  of  the  material  Ftu  is  used.  This  applied  load  Pappiied  is  given  by 

P applied  =  Cu  D  ts 

where  au  is  either  the  design  ultimate  tensile  stress  or  the  tensile  ultimate  skin 
strength,  ts  is  the  skin  thickness,  D  is  the  diameter  of  antenna  connector  cutout 
hole.  The  margin  of  safety  is  given  by 


P 

Margin  of  Safety = — ssi_-1 

Pappiied 

2.3  The  Shear  Margin  of  Safety 

The  shear  margin  of  safety  of  the  antenna  installation  is  calculated  by 


Margin  of  Safety  =  - 1 


where  Fsu  is  the  ultimate  shear  strength  of  the  skin  or  mounting  plate  material. 
The  equation  is  used  for  each  side  of  the  mounting  plate. 

2.4  Margin  of  Safety  as  a  Criterion 

The  margins  of  safety  (MS)  based  on  the  mounting  plate  allowable  and  the 
fastener  joint  allowable  were  calculated  to  determine  the  adequacy  of  the 
antenna  installation  on  fuselage  skin. 

•  MS  <  0  Antenna  installation  is  statically  inadequate 

•  MS  =  0  Antenna  installation  is  marginally  adequate 

•  MS  >  0  Antenna  installation  is  statically  adequate 

For  antenna  installations  that  are  not  statically  adequate,  they  must  be 
redesigned  to  ensure  the  adequacy  of  the  antenna  installation. 
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2.5  The  Stiffness  Check  of  the  Antenna  installation 

The  stiffness  ratio  between  the  mounting  plate  and  the  skin  is  calculated  using 
the  following  equation 


Stiffness  Ratio  = - - 

C Et)s 

The  antenna  installation  is  considered  adequate  if  the  ratio  is  between  1 .0  and 
1 .5.  The  antenna  installation  is  too  stiff  when  the  value  is  greater  than  1 .5  and 
not  stiff  enough  when  it  is  less  than  1 .0. 

2.6  The  Fastener  Bending  Check  of  the  Antenna  Installation 

The  fastener  bending  is  checked  using  the  following  equation 


d 

Where  d  is  the  fastener  diameter,  ts  and  tp  are  the  thickness  of  skin  and  mounting 
plate,  respectively.  The  parameter  Q  is  the  fastener  bending  indicator.  For 
aluminum  fasteners,  the  bending  is  important.  A  Q  value  above  2  may  indicate 
that  the  aluminum  rivet  will  not  fill  the  hole  but  instead  may  buckle  in  the  hole.  In 
such  a  case,  RAPIDC  recommends  steel  or  titanium  fasteners  be  used.  For 
steel  and  titanium  fasteners,  there  is  no  constraint  for  typical  fuselage  antenna 
installations. 

2.7  The  Inter-Rivet  Buckling  Guideline 

To  avoid  inter-rivet  buckling  in  the  modified  skin,  the  fastener  spacing  should  be 
four  to  six  times  the  diameter  of  the  fastener  shank  diameter.  This  guideline  is 
enforced  during  user  input  phase  of  RAPIDC. 

3.0  Damage  Tolerance  Analysis 

To  perform  the  damage  tolerance  analysis  of  a  modified  skin  for  antenna 
installations,  critical  locations  in  the  skin  must  first  be  determined.  Assumptions 
of  initial  flaws  at  fracture  critical  locations  and  the  continuing  damage  need  to  be 
made.  The  stress  spectrum  must  also  be  prescribed.  In  addition,  the  following 
data  are  needed: 

•  Crack  growth  rate  data  of  the  skin  material 

•  Fracture  toughness  of  the  skin  material 

•  Stress  intensity  factors  of  relevant  crack  configurations 
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Damage  tolerance  analysis  procedure  is  described  in  the  flowchart  shown  in  - 
Figure  2.  Each  element  in  the  analysis  procedure  is  described  as  follows. 


No 


Critical  Fastener 
Location  &  Load 


Initial  Crack 
Geometry 


± 


Stress  Intensity 
Factors 


Inspection 
Threshold  &  Interval 


Crack 

Advance 

Crack 

Growth 

J 

i 

Yes 

_ V 

Residual 

— ^  j 

Strength 

Stress 

Spectrum 


Crack  Growth 
Rate  Data 


Figure  2.  Flow  Chart  of  Damage  Tolerance  Analysis  Procedure 


3.1  Critical  Fastener  Locations 

Damage  tolerance  analysis  of  the  modified  skin  begins  by  postulating  the  initial 
flaws  at  the  critical  locations.  For  rectangular  installations,  two  locations  are 
identified  at  the  center  and  corner  fastener  holes  in  the  outermost  fastener  rows 
that  have  higher  load  transfer.  A  third  location  is  identified  at  the  edge  of 
antenna  connector  hole  cutout  in  the  skin  where  stress  concentration  is  the 
highest.  For  circular  installations,  the  critical  locations  are  identified  at  the  top  or 
bottom  fastener  in  the  outermost  ring  and  at  the  antenna  connector  hole  cutout 
as  well.  The  critical  fastener  locations  in  circular  and  rectangular  installations  are 
indicated  in  Figure  3. 
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Figure  3.  Critical  Fastener  Locations 


3.2  Fastener  Loads  and  Skin  Stress  Gradients 

The  fastener  loads  or  skin  stress  gradients  along  the  potential  crack  growth  path 
are  needed  for  calculating  crack  tip  stress  intensity  factors.  The  determination  of 
fastener  loads  and  skin  stress  gradients  was  accomplished  using  the 
FRANC2D/L  program  developed  jointly  at  Cornell  University  and  Kansas  State 
University  sponsored  by  NASA-Langley  Research  center.  The  finite  element 
model  was  created  using  the  FRANC2D/L  preprocessor  program  called  CASCA 
developed  at  Cornell  University. 

In  FRANC2D/L  analysis,  the  fastener  shear  rigidity  needs  to  be  calculated 
separately  and  input  into  the  program.  The  FRANC2D/L  program  then 
assembles  stiffnesses  contributing  from  skin,  mounting  plate,  and  fasteners  to 
form  the  overall  stiffness  matrix.  The  shear  rigidity  can  be  calculated  using  the 
Swift’s  empirical  equation: 


K= _ Ed _ 

A+B(— + - - - ) 

t  t 

lSkin  1  Mounting  PI  ate 

where  E  is  the  weighted  Young’s  modules  of  the  skin  and  mounting  plate,  d  is 
the  fastener  hole  diameter,  tSkin  and  tMounting  Plate  are  the  thickness  of  the  skin  and 
mounting  plate,  respectively,  and  A  and  B  are  material  dependent  empirical 
constants.  The  values  of  A  and  B,  for  aluminum  fasteners,  are  5.0  and  0.8, 
respectively. 

In  RAPIDC,  databases  of  fastener  loads  and  skin  stress  gradients  along  the 
potential  crack  growth  path  are  developed  for  circular  and  rectangular 
installations.  Parameters  considered  include: 
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Skin  thickness:  0.032”,  0.040”,  0.050",  and  0.063” 

Mounting  plate  thickness:  1  gauge  higher  than  skin  thickness 
Cutout  diameter:  1 .0”,  2.0”,  and  3.0" 

Fastener  size:  1/8”,  5/32”,  and  3/16” 

Number  of  fastener  rows:  1  ~  5 
Aspect  ratio:  1 :1  (square),  1 :1 .3  and  1 :1 .6 
Materials:  Aluminum  skin  and  mounting  plate 
Aluminum  fastener 

Fastener  loads  and  skin  stress  gradients  calculations  along  the  potential  crack 
growth  paths  are  described  in  Appendix  A. 

3.3  Initial  Flaw  Size  and  Subsequent  Growth  Scenarios 

To  simplify  the  damage  tolerance  analysis  of  the  modified  skin;  initial  flaws  and 
their  growth  pattern  were  assumed.  Three  scenarios  are  considered  in  RAPIDC 
for  rectangular  antenna  installations,  according  to  the  decision  made  in  a 
technical  review  meeting  (September  22,  1999  in  Albuquerque,  New  Mexico), 
and  are  described  below  in  Figure  4. 

Scenario  1 :  Center  fastener  hole  in  the  outermost  fastener  row 

Initial  Crack:  Two  diametric  through  cracks  of  lengths  0.05”  and  0.005”, 
respectively,  emanating  from  the  center  fastener  hole 
together  with  a  0.005”  crack  at  one  side  of  every  other  hole 


O  o0  005^)0005^:05" 


0.005"  0.005” 

c>  o  o 


Subsequent  Damage:  All  cracks  grow  concurrently  but  independently, 
interaction  between  cracks  being  ignored.  The  amount  of 
growth  8ai  for  the  0.005”  crack  is  added  to  its  original  length 
when  the  0.05"  crack  grows  into  the  adjacent  hole.  The 
same  process  continues  in  successive  growth. 

0.005”+  Sa1  0.005”+  5a, 

~o  ~o  -o  -O — o~  o~  o 

0.005”+  5a,  +  5a2  0.005”+  5a,  +  Sa2 

■O  -O  "O — O — O — O-  O- 
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Scenario  2:  Corner  fastener  hole  in  the  outermost  fastener  row 

Initial  Crack:  Two  diametric  through  cracks  of  lengths  0.05”,  pointing 
toward  the  adjacent  hole,  and  0.005”,  respectively, 
emanating  from  the  corner  fastener  hole  together  with  a 
0.005”  crack  at  one  side  of  every  other  hole 


O  O  O 


0.005"  0.005”  0.05”  0.005” 


Subsequent  Damage:  All  cracks  grow  concurrently  but  independently, 
without  interaction  between  cracks  being  considered.  The 
amount  of  growth  5a  for  the  0.005”  crack  is  added  to  its 
original  length  when  the  0.05”  crack  grows  into  the  adjacent 
hole.  The  same  process  continues  in  successive  growth. 

0.005”+  5a1L  0.005”+  5alR 

O  -O  -O  -O  -O — O 

0.005”+  5cIil  +  Sa2i  0.005”+  53ir  +  Sa2  r 

0  0  0-0  0—0 — o- 


It  is  noted  that  successive  growths  of  the  0.005”  at  corner 
fastener  hole  Sam,  5a2R,  etc.  are  different  from  the  growths  at 
other  holes  5am  6a2i_, etc.  because  the  former  does  not  grow 
toward  an  adjacent  hole  but  others  do. 

Scenario  3:  Antenna  connector  hole 

Initial  Crack:  Two  diametric  through  cracks  of  lengths  0.05”  and  0.005”, 
respectively,  emanating  from  the  antenna  connector  hole 
together  with  a  0.005”  crack  at  one  side  of  every  other  hole 


0.005”  0.005”  0.05"  ' 0.005”  0.005” 

O  O  \_J-  o-oo- 


Subsequent  Damage:  All  cracks  grow  concurrently  but  independently, 
ignoring  interaction  between  cracks.  The  amount  of  growth 
•5a  for  the  0.005”  crack  is  added  to  its  original  length  when 
the  0.05”  crack  grows  into  the  adjacent  hole.  The  same 
process  continues  in  successive  growth. 
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It  is  noted  that  successive  growths  of  the  0.005”  at  antenna 
connector  hole  San.,  5a2L,  etc.  are  different  from  the  growths  of 
0.005”  crack  at  fastener  holes  Sais.  Sa2s,  etc.  because  the  hole 
sizes  and  fastener  pitches  are  different. 

Figure  4.  Initial  Flaws  and  Subsequent  Growths  Assumptions 

For  circular  antenna  installations,  two  scenarios  are  considered.  The  first 
scenario  postulates  two  initial  diametric  cracks  with  0.05”  and  0.005”  at  the  top  or 
bottom  fastener  in  the  outermost  ring.  The  cracks  grow  concurrently  but 
independently  in  the  skin.  The  second  scenario  is  the  same  as  Scenario  3  for 
the  rectangular  antenna  installation. 

3.4  Stress  Intensity  Factors 

In  RAPIDC,  stress  intensity  factors  are  calculated  using  three  methods.  For  a 
crack  emanating  from  a  pin-loaded  fastener  hole  in  the  skin,  the  stress  intensity 
factor  is  calculated  from  the  gross,  bearing,  and  bypass  stresses  using  the 
superposition  method.  For  a  crack  emanating  from  the  antenna  connector  open 
hole  in  the  skin,  it  is  calculated  from  stress  gradients  using  the  weight  function 
method.  To  account  for  the  hole  effect  as  a  crack  is  growing  toward  the  adjacent 
hole,  an  engineering  approach  which  compounds  geometry  factors  is  used. 
These  methods  are  described  below. 

Method  1 :  Superposition  method 

As  the  method  implied,  the  stress  intensity  factor  for  a  crack  emanating  from  a 
hole  in  a  plate  subjected  to  gross,  bearing,  and  bypass  stresses  can  be  obtained 
using  the  superposition  method  as  illustrated  in  Figure  5. 
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**  Gross 

t  . 

r 

^  Gross  *  ®  Bypass 

Gki"1 

^  Bearing 

=  1/2  < 

(“on 

+ 

f  ®  Bearing 

(Si 

^  &  Bearing 

'■v 

> 

^  Bypass 

1 

a  Gross  +  &  Bypass 

J 

Figure  5.  Superposition  of  Stress  Intensity  Factors 

Let  BRF  and  BPF  denote  the  bearing  factor  (asearing  /  across  )  ancl  the  bypass 
factor,  (cfBypass  /  across  ),  respectively.  The  stress  intensity  factor  K  can  be 
calculated  by 


K  K  Far  Field  K  Point  Loads 

=  A 0+  BPF)aGross Vira  PFafField  +  X  BRF  aGrossVrca  PPointLoads 

=  X  [(l  +  BPF)pFarField  +  BRF  pPoimLoads  ]aGross  4k a 

In  the  above  equation,  pFar  Reid  and  pPOint  Loads  are  the  geometry  factors  (or 
P  factors)  which  are  the  normalized  stress  intensity  factors  of  a  crack  at  a  hole  in 
a  wide  plate  subjected  to  a  unit  far  field  stress  and  unit  bearing  stress  divided  by 
4Ha ,  respectively.  The  crack  length  “a”  is  measured  from  the  crack  tip  to  the 
edge  of  the  hole. 

The  stress  intensity  factor  K  is  divided  by  Koto  obtain  the  geometry  factor  p  as 

P  -  K/Ko 


=  Va  [  (1  +  BPF  )  pFar  Field  +  BRF  pPoint  Loads  ] 

where  Ko  is  the  stress  intensity  factor  of  a  crack  at  a  hole  in  a  wide  plate 
subjected  to  a  gross  stress  across,  i.e.,  K0  =  aGross  4ksl . 

To  calculate  the  geometry  factor  p,  the  geometry  factors  pFar  Reid  and  pPOim  Loads 
are  needed.  Shown  in  Figures  6  and  7  are  geometry  factors  pFar  Reid  and  pPoint 
Loads,  respectively. 
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wWiiiFi 


Jized  Crack  Length:  a/R 


re  6.  A  Crack  Emanating  from  a  Hole  under 
Uniform  Far  Field  Stress:  K0  =  a 


Figure  7.  A  Crack  Emanating  from  a  Hole  under  a  Pair 

F  _ _ 

of  Point  Loads:  K0  =  — —  Vrca 
0  2Rt 


Method  2:  Weight  function  method 

The  stress  intensity  factor  of  a  crack  emanating  from  a  hole  in  a  wide  plate 
subjected  to  stress  gradients  on  crack  faces  as  shown  in  Figure  8(a)  can  be 
calculated  using  the  weight  function  method.  In  this  case,  the  stress  intensity 
factor  is  obtained  by  summing  up  stress  intensity  factors  for  each  pair  of  point 
loads  acting  on  crack  faces  at  a  distance  x  from  the  edge  of  the  hole  as  shown  in 
Figure  8(b). 


(a)  Stress  Gradients 


(b)  A  Pair  of  Point  Loads 


Figure  8.  Illustration  of  the  Weight  Function  Method 


Suppose  the  stress  intensity  factor  of  the  crack  geometry  subjected  to  the  load 
condition  described  in  Figure  8(b)  is  known  as 

AK(p,a,x)  =  pyfOa  G(x) 

where  G(x)  is  the  Green’s  function  representing  the  geometry  factor  of  the  crack 
geometry  under  a  unit  of  point  loads  on  crack  faces.  The  stress  intensity  factor 
solution  for  the  problem  in  Figure  8(a)  can  be  calculated  as 

K(p,a)  =  J^AK(p,a,x)dx 

For  the  stress  gradient  cr(x)  =  cro  <j>(x)  as  shown  in  Figure  8(a),  the  load  per  unit 
thickness  at  a  distance  x  from  the  edge  of  the  hole  can  be  expressed  by 

p(x)  =  cy0<|>(x)dx 

Therefore,  the  stress  intensity  factor  for  the  same  geometry  subjected  to  the 
stress  gradient  cr(x)  =  ct0<Kx)  on  crack  faces  can  be  calculated  as 

K  =  cr0  -Jna.  <J>(x)G(x)dx 

The  geometry  factor  p  can  be  obtained  as 

P  =  K/K0 

where  K0  is  the  stress  intensity  factor,  i.e.,  K0  =a0 Tjca . 
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Method  3:  Compounding  method 


To  account  for  the  hole  effect  on  the  stress  intensity  factor  of  a  crack  at  a  hole 
growing  toward  an  adjacent  hole,  an  engineering  approach  by  compounding  the 
relevant  geometry  factors  is  employed.  The  method  can  be  illustrated  using  an 
example. 

Suppose  the  crack  tip  stress  intensity  factor  of  the  geometry  shown  in  Figure  9  is 
sought. 


Figure  9.  A  Crack  at  a  Hole  Approaching  an  Adjacent  Hole 


The  stress  intensity  factor  can  be  estimated  as 

K  —  C7q  -Jk  a  P  Compound 

where  pcomtpound  is  the  compounded  geometry  factor  calculated  as  the  product  of 
the  geometry  factors  (3At  a  Hole  and  PApproaching  a  Hotel  i*®*» 

Pcompounded  =  PAt  a  Hole  X  PApproaching  a  Hole 

For  the  case  of  a  crack  emanating  from  a  hole  in  a  wide  plate  subjected  to  stress 
gradients  on  crack  faces  as  shown  in  Figure  1 0,  the  geometry  factor  pAt  a  Hole  can 
be  obtained  using  the  weight  function  method  described  in  Method  2.  For  the 
case  of  a  crack  growing  toward  a  hole,  the  geometry  factors  pAPProaching  a  Hole  are 
obtained  from  an  open  literature  and  are  presented  in  Figure  11. 
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Figure  10.  A  Crack  Approaching  a  Hole 
K0  =0*7*1 


0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 

Normalized  Crack  Length:  a/b 


Figure  1 1 .  Geometry  Factors  for  a  Crack  Approaching  a  Hole 

3.5  Residual  Strength  Evaluation  and  Limit  Stress 
3.5.1  Fracture  Toughness  Method 

The  residual  strength  of  the  modified  skin  for  antenna  installation  is  calculated 
using  the  fracture  toughness  method 
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where  Kc  is  the  toughness  of  the  skin  material,  p  is  the  geometry  factor,  and  a 
denotes  the  crack  length.  For  small  crack  lengths,  the  residual  strength  is 
determined  using  the  Fedderson’s  method.  In  this  method,  a  tangent  line  from 
the  yield  stress  ay  of  the  material  to  the  residual  strength  curve  obtained  using 
the  fracture  toughness  method  is  numerically  obtained  as  illustrated  below  in 
Figure  12. 


Figure  1 2.  Residual  Strength  Curve 


3.5.2  Limit  Stresses 


To  assess  the  residual  strength  capability  of  the  modified  skin,  the  residual 
strength  calculated  in  Section  3.5.1  is  checked  against  the  limit  stress  of  the  skin. 
In  the  analysis,  the  crack  continues  to  grow  until  the  residual  strength  of  the  skin 
is  no  longer  greater  than  the  limit  stress. 


The  limit  stress  in  the  circumferential  direction,  for  the  longitudinal  crack,  is 
calculated  using  the  equation 


'  Limit, Circumferential 


=  1.1  — +  0.6  — 


t 


t 


where  p  is  the  pressure  differential,  R  and  t  are  the  radius  and  thickness  of  the 
fuselage  shell  at  the  location  where  the  damage  tolerance  analysis  is  performed. 
The  stress  due  to  the  additional  pressure  0.6  psi  in  the  equation  is  to  account  for 
the  external  aerodynamic  suction  pressure  (0.5  psi)  plus  a  0.1  psi  regulator 
tolerance  per  FAR  23.574. 

The  limit  stress  in  the  longitudinal  direction,  for  the  circumferential  crack,  is 
determined  using  the  equation 
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pR  ..R  , 

^Limitlongitudinal  ~  ^  “h  ^  ^ 

where  p  is  the  pressure  differential,  n,  is  the  maneuvering  limit  load  factor,  y  is 
the  payload  reduction  factor,  R  and  t  are  the  radius  and  thickness  of  the  fuselage 
shell,  and  ctig  is  the  one-G  stress  at  the  location  where  the  damage  tolerance 
analysis  is  performed.  The  additional  pressure  0.6  psi  is  added  per  FAR  23.574. 


3.6  Load/Stress  Spectrum  and  Equivalent  Stress 

To  perform  the  cycle-by-cycle  crack  growth  analysis  of  a  modified  skin  to 
antenna  installation,  the  stress  spectrum  that  exerts  on  the  skin  needs  to  be 
prescribed.  The  stress  spectrum  is  generally  derived  from  the  airplane  C.G.  load 
spectrum  which  is  developed  from  the  airplane  flight  profiles,  exceedance  data, 
and  pressure  schedule.  The  procedures  of  load  and  stress  spectra  development 
for  commuter  airplanes  are  described  in  Appendix  B. 

In  RAPIDC,  the  stress  spectrum  can  be  converted  into  a  constant  amplitude 
stress  spectrum  representing  a  single  cycle  of  an  equivalent  repeated  flight.  The 
constant  amplitude  stress  spectrum  is  then  used  in  the  simplified  crack  growth 
analysis.  The  procedure  is  described  in  Appendix  C. 

3.7  Crack  Growth  Rate  Data 

The  crack  growth  da/dN  vs.  AK  data  are  obtained  experimentally  from  constant 
amplitude  coupon  tests  for  the  material,  and  are  normally  documented  in  a 
tabular  or  graphic  form.  The  tabulated  data  can  be  used  directly  in  the  crack 
growth  calculation  using  the  cycle-by-cycle  method.  To  use  the  simplified 
method,  the  tabular  data  need  to  be  curve-fitted  to  the  Walker's  equation: 

fj-c[(l-R)'K_f 

In  RAPIDC,  the  following  da/dN  versus  AK  tabular  data  as  well  as  the 
coefficients  C,  p,  and  q  in  the  Walker’s  equation  of  thirteen  materials  are 
included  in  the  material  database. 


(1 )  2024-T3  Clad  Sheet,  -T42  Bare  Sheet,  L-T  RT  LA  DW 

(2)  2024-T3  Clad  Sheet,  -T42  Bare  Sheet,  T-L  RT  LA  DW 

(3)  2024-T351  Plate,  -T351 1  Extrusion,  L-T  RT  LA  DW 

(4)  7050-T7452  Forging,  L-T  T-L  LA  RT 

(5)  7050-T7451 1 ,  -T7651 1  Extrusion,  L-T  RT  LA  DW 

(6)  7050-T7451 1  Extrusion,  L-T  RT  STW 

(7)  7050-T765 1 ,  -T745 1  Plate,  L-T  T-L  RT  LA  DW 
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(8)  7050-T7651 1  Extrusion,  L-T  RT  STW 

(9)  7475-T7351  Plate,  L-T  LA  RT 

(1 0)  7475-T7651  Plate,  L-T  LA  RT  DW 

(1 1 )  7475-T761  Clad  Sheet,  L-T  RT  LA  DW 

(1 2)  7075-T6  Clad  Sheet,  L-T  RT  LA 

(1 3)  201 4-T6  Sheet,  T=0.05  L-T  RT  LA  1 0  HZ 

where 

L-T :  longitudinal-transverse  material  orientation 
T-L:  transverse-longitudinal  material  orientation 
RT:  room  temperature  test  condition 
LA:  laboratory  air  environmental  test  condition 
DW:  distilled  water  environmental  test  condition 
STW:  sump  tank  water  environmental  test  condition 

RAPIDC  also  provides  user  with  a  direct  input  of  da/dN  tabular  data  or  Walker’s 
coefficients.  The  input  format  of  the  da/dN  tabular  data  is  described  as  follows. 


Line  1  1 


3  30  (2024-T3  CLAD  SHEET, -T42  BARE  SHEET,  L-T  RT  LA  DW) 

0.05  0.40  0.80 


Line  4*>O.IOOE-3I 

0.100E-07 

0.300E-07 

0.500E-07 

0.700E-07 

0.800E-07  0.100E-06 

0.200E-06 

0.300E-06 

0.500E-06 

0.800E-06 

0.100E-05 

0.200E-05 

0.500E-05  0.800E-05 

0. 1 00E-04 

0.200E-04 

0.400E-04 

0.800E-04 

0.100E-03 

0.200E-03 

0.500E-03  0.100E-02 

0.300E-02 

0.800E-02 

0.200E-01 

0.500E-01 

0.100E+00  0.200E+00  0.100E+01 

Line  ->18990 

2.9000 

2.9010 

3.0000 

3.2000 

3.3000  3.6000 

4.7500 

5.5000 

6.2000 

6.6000 

6.8000 

8.0000 

10.1000  12.2000 

13.2500 

16.9000 

20.0000 

24.2000 

25.6000 

29.6000 

36.0000  41.0000 

50.5000 

61.0000 

70.0000 

78.0000 

87.0000 

92.0000 

93.0000 

Line  1 2>2.3480 

2.3490 

2.3500 

13600 

2.5000 

2.6000  2.8800 

3.9000 

4.5000 

5.1500 

5.7000 

6.0000 

7.0000 

8.7500  9.9000 

10.5000 

13.2500 

16.0000 

19.5000 

21.0000 

24.5000 

29.0000  33.5000 

41.5000 

48.0000 

51.0000 

54.0000 

56.5000 

57.0000 

58.0000 

Line  1&>i.7490 

1.7500 

1.7600 

1.7700 

1.8200 

1.8750  2.0500 

2.8500 

3.4000 

4.0000 

4.5000 

4.8000 

5.6000 

7.0000  8.0000 

8.4000 

9.9000 

11.5000 

13.0000 

13.5000 

15.3000 

18.0000  19.8000 

20.0000 

20.0010 

20.0020 

20.0030 

20.0040 

20.0050 

20.0060 
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Definition 


Line 

1 


2 


2 


2 

3 


4-7 


8-11, 

12-15, 

16-19 


Acronym 

NMAT 


MRATIO 


MVALUE 


RATIOV 


DADN 


DELTAK 


Type 

Integer 


Integer 


Integer 


Real 


Real 


Real 


The  material  number  that  corresponds  with  the 
material  number  containing  the  Walker’s  C,  p 

and  q  values _ _ _ _ 

The  first  integer  number  is  the  number  of 
different  ratio  values  in  line  3;  the  maximum 

value  is  10 _ _ _ 

The  second  integer  value  is  the  number  of  da/dN 

values  in  lines  4-7,  and  AK  values  in  lines  8-1 1 , 

12-15,  and  16-19;  the  maximum  value  is  40 

Description  of  the  material;  optional 

The  different  ratio  values  in  ascending  order; 

there  are  MRA 770  values 

The  da/dN  values;  there  are  MVALUE  values  for 

this  v-axis  _ _ 

The  AK  values;  there  are  MVALUE  values,  and 
MRATIO  sets  for  this  x-axis 


3.8  Crack  Growth  Analysis 

Two  methods,  the  simplified  method  and  the  cycle-by-cycle  method,  are 
implemented  in  RAPIDC  to  calculate  the  crack  growth  life  of  the  modified  skin. 

3.8.1  Simplified  Method 

The  crack  growth  of  a  cracked  skin  can  be  calculated  using  the  simplified 
method  based  on  the  Walker's  crack  growth  equation.  The  number  of  flights  Njj 
for  a  crack  growing  from  the  size  to  the  size  a.  (a.  >  aj)  is  calculated  as 

'V£<S°ij>"P 


where  S  is  the  one  cycle  equivalent  stress,  C  and  p  are  the  coefficients  in 
Walker's  equation,  and  the  geometry  term  Gij'p  can  be  calculated  as 


GjfP  =  J3‘  [p(a)V^a]  pda 


Let  the  integrand  in  the  above  integral  be  represented  by  the  geometric  function 
g(a),  i.e., 

g(a)  =  [|3(a)  Vrea"  ]"p 


where  P(a)  is  the  compounded  normalized  stress  intensity  factor  obtained  as 


P(a)  P  AtorPassingTtiroughHole(s)  P  ApproachingHole 


The  terms  PAt  or  Passing  Through  Hole(s)  Slid  ^Approaching  Hole  STe  the  geometry  factors  Of 

relevant  crack  configurations  as  described  in  Sections  3.4. 

The  integral  of  the  geometric  function  represents  the  geometry  terms.  Since  the 
geometry  factors  are  normally  presented  in  a  tabular  form,  the  integral  is  carried 
out  numerically  using  the  Gaussian  Quadrature  method. 

In  the  Gaussian  Quadrature  method,  the  integral  can  be  approximated  by  the 
following  equation: 


where 


Gii'p=“1j^XWk9(ak) 

k=1 

g(ak)  =  [p(ak)V^rp 


and 


ak 


(aj  -  a,)tk  +  (a,-  +  a;) 
2 


In  the  above  equations,  the  variables  wk  and  tk  are  the  k-th  weighting  coefficient 
and  root  of  the  K-th  order  Legendre  polynomial  P«  (t)  =  0.  For  instance, 
numerical  values  of  w  and  t  for  the  case  of  K  =  4  are  given  below. 


k 

wk 

tk 

1 

0.347855 

-0.861136 

2 

0.652145 

-0.339981 

3 

0.652145 

0.339981 

4 

0.347855 

0.861136 

It  is  noted  that  the  simplified  method  calculates  the  number  of  flights  Nq  for  a 
crack  growing  from  a  size  a.  to  a  next  size  aj  (aj  >  aj).  With  this  method,  the 
number  of  flights  can  be  calculated  in  one  step  for  a  crack  growing  from  an  initial 
size  to  the  final  size.  To  generate  crack  growth  life  at  intermediate  crack  lengths, 
the  procedure  is  carried  out  between  two  crack  lengths  and  iterated  to  complete 
the  crack  growth  analysis. 
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3.8.2  Cycle-by-Cycle  Method 

The  cvcle-by-cyc!e  crack  growth  analysis  method  is  illustrated  using  an  example 
as  shown  in  figured  A  flow  chart  of  the  analysis  procedure  is  provided  in 

Figure  14. 


a 

(a)  Crack  at  hole  in  a  (t>)  Stress  Spectrum 


Normalized  Crack  Length:  a/R  °  Stress  Intensity  Factor  aK 


Figure  13.  Description  of  Crack  Growth  Analysis  Method 
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Figure  14.  Flow  Chart  of  the  Cycle-by-Cycle  Crack  Growth  Analysis 


Let  ao  be  the  initial  crack  length.  The  crack  growth  analysis  starts  with  an  initial 
crack  of  length  ao  at  i  =  1 .  For  the  first  cycle  stress  excursion,  the  stress-intensity 
factor  and  the  plastic  zone  size  are  calculated.  With  the  stress  ratio  R  and  the 


change  of  the  stress-intensity  factor  AK,  the  incremental  crack  length  is  obtained ' 
from  the  da/dN  data.  The  crack  is  incremented  to  the  length  ai. 

At  each  sequential  i-th  cycle,  the  crack  growth  starts  with  a  length  equal  to  an. 
The  stress-intensity  factor  and  the  plastic  zone  size  are  calculated  using  the  i-th 
stress  excursion.  A  test  is  then  made  to  determine  if  the  growth  of  the 
progressing  crack  is  retarded.  A  positive  answer  leads  the  analysis  to  use  the 
Generalized  Willenborg  retardation  model  in  the  calculation  of  incremental  crack 
length.  Otherwise,  the  retardation  effect  is  bypassed  in  the  calculation.  During 
the  i-th  cycle,  the  crack  advances  a  length  of  &a{  and  grows  to  a  length  a\  =  au  + 
Aaj. 

The  process  continues  until  the  crack  grows  to  the  critical  crack  length 
determined  from  the  residual  strength  analysis. 

In  the  analysis,  the  compression-tension  stress  cycle  is  treated  as  a  zero-tension 
stress  cycle  and  the  compression-compression  stress  cycle  is  ignored  in  the 
crack  growth  calculation.  Furthermore,  the  parameter  <)>  introduced  in  the 
generalized  Willenborg  retardation  model  to  calculate  the  effective  stress 
intensity  factor  (see  flowchart  in  Figure  14)  is  calculated  by  the  equation 

■j  _  AK  Threshold 


where  AKihreshoid  is  the  threshold  stress-intensity  factor  level  associated  with  zero 
fatigue  crack  growth  rates,  and  S0L  is  the  overload  (shut-off)  ratio  required  to 
cause  crack  arrest  for  the  given  material.  A  value  of  2.3  is  used  in  RAPIDC. 

4.0  Inspection  Threshold  and  Interval 

With  the  crack  growth  life  of  the  modified  skin  provided  by  RAPIDC,  the  user 
determines  the  inspection  threshold  and  the  inspection  interval  at  the  detectable 
crack  length  for  the  desired  inspection  method. 


5.0  Example  Problem 

The  damage  tolerance  analysis  of  a  modified  skin  to  an  antenna  installation  is 
performed  below  for  three  crack  scenarios  described  in  Section  3.3. 

Description  of  the  Example  Problem 

Figure  15  depicts  an  example  problem  of  an  antenna  installation  on  the  fuselage 
of  a  commuter  airplane.  Assume  that  the  skin  is  subjected  to  an  equivalent  one 
cycle  circumferential  stress  of  18  ksi  with  zero  stress  ratio.  Damage  tolerance 
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analysis  of  the  modified  skin  is  carried  out  to  demonstrate  the  analysis - 
procedure. 


Fuselage  radius:  60.0" 

Skin  thickness:  0.050” 

Mounting  plate  thickness:  0.063” 
Cutout  diameter:  1 .0” 

Fastener  diameter:  3/16" 

Fastener  rows:  3 
Fastener  pitch:  1.0" 

Skin  and  mounting  plate:  AI2024T3 
Fastener:  NASI  097 
Fracture  toughness:  153  ksi  sqrt(in) 
Pressure  differential:  8.6  psi 


Figure  15.  Description  of  the  Example  Problem 

Fastener  Loads  and  Stress  Gradients 

(1)  Fastener  Loads 

Figure  16  shows  load  transfers  at  fasteners  in  pounds  due  to  1000  psi 
reference  skin  stress.  These  fastener  loads  were  obtained  using  FRANC2DL 
as  described  in  Appendix  A.  As  expected,  the  outer  row  is  the  critical 
fastener  row. 
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Figure  16.  Fastener  Loads 


2.  Stress  Gradients 

The  crack  growth  scenario  3  described  in  Section  3.3  is  depicted  in  Figure  17. 
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Potential  Crack 
Growth  Path. 
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Figure  1 7.  Potential  Crack  Growth  Paths 


The  stress  gradient  along  the  potential  crack  path  subjected  to  1,000  psi 
reference  skin  stress  is  shown  in  Figure  18.  The  distance  d  is  measured  from 
the  center  of  the  hole. 


Figure  18.  Stress  Gradient  along  Potential  Crack  Growth  Path 


Residual  Strength  and  Crack  Growth 

To  demonstrate  the  residual  strength  and  crack  growth  analysis  of  the  example 
problem,  analysis  is  performed  for  three  Scenarios  described  in  Section  3.3. 
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(1)  Stress  Intensity  Factors 

Stress  intensity  factors  for  a  crack  initiating  from  the  center  or  the  corner  fastener 
hole  due  to  gross,  bearing,  and  bypass  stresses  are  calculated  using  the 
superposition  method.  The  effect  of  a  crack  growing  toward  an  adjacent  hole  is 
accounted  for  using  the  compounding  method. 

For  a  crack  initiating  from  the  antenna  connector  hole,  the  weight  function 
method  is  used.  Using  the  stress  gradients  shown  in  Figure  18,  the  normalized 
stress  intensity  factor  or  geometry  factor  is  obtained  as  shown  in  Figure  1 9. 


Distance  d  (in) 


Figure  19.  Normalized  Stress  Intensity  Factors 

The  effect  of  a  crack  growing  toward  an  adjacent  hole  is  accounted  for  using  the 
compounding  method. 
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(2)  Residual  Strength 


Residual  strengths  of  the  skin  are  calculated  using  the  fracture  toughness 
method  and  modified  using  Fedderson’s  criterion.  The  calculation  terminates 
when  the  residual  strength  reaches  the  limit  stress  of  12,072  psi  obtained  using 
the  circumferential  limit  stress  equation  in  Section  3.5.2.  Figure  20  shows  the 
residual  strength  curve. 


Figure  20.  Residual  Strength 


(3)  Crack  Growth  Lives 

Figure  21  shows  the  crack  growth  lives  obtained  using  the  simplified  crack 
growth  analysis  method. 


Figure  21 .  Crack  Growth  Life 
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Inspection  Threshold  and  Interval 


Once  the  residual  strength  and  crack  growth  life  have  been  obtained,  the  critical 
crack  length  can  be  determined  from  the  residual  strength  at  the  limit  stress. 
The  crack  growth  life  at  the  critical  crack  length  can  then  be  determined  from  the 
crack  growth  life  curve. 

With  the  crack  growth  history,  the  inspection  threshold  and  interval  can  be 
determined  from  the  crack  growth  life  based  on  the  detectable  crack  size  of  the 
inspection  method. 
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Appendix  A 

Fastener  Loads  and  Skin  Stress  Gradients  Calculations 
A1.  Introduction 

This  appendix  describes  the  database  development  of  load  transfers  along  the 
critical  outermost  fastener  row  and  stress  gradients  along  the  potential  crack 
growth  path,  for  a  flaw  initiating  from  the  antenna  connector  hole,  normal  to  the 
stress  application  direction  in  a  modified  skin  to  antenna  installation.  The 
purpose  of  the  work  is  to  determine  fastener  loads  and  stress  gradients  required 
in  the  calculation  of  stress  intensity  factor  for  the  damage  tolerance  analysis. 
The  database  is  developed  for  square/rectangular  and  circular  antenna 
installations  typically  found  in  commuter  airplanes.  A  similar  database  will  be 
developed  for  elliptical  installations  in  the  next  phase  of  contract. 

A2.  Types  of  Antenna  Installation 

Three  types  of  antenna  installation  are  common  on  the  fuselage  skin  of 
commuter  airplanes.  Typical  mounting  plates  in  most  antenna  installations 
including  rectangular,  circular,  and  elliptical  shapes  are  shown  in  Figure  A1 .  In 
all  installations,  the  mounting  plate  is  mechanically  fastened  to  the  skin  around 
the  antenna  connector  hole  cutout. 


(a)  Rectangular 


(b)  Circular 


(c)  Elliptical 


Figure  A1 .  Three  types  of  antenna  installation 

In  RAPIDC  Version  1 .0,  only  the  circular  and  square/rectangular  installations  are 
considered.  The  elliptical  installations  will  be  included  in  the  next  phase  of 
contract. 


A3.  Approach 

To  perform  the  damage  tolerance  analysis,  fastener  loads  or  skin  stress 
gradients  along  the.  potential  crack  growth  path  are  needed  when  calculating 
crack  tip  stress  intensity  factors.  The  determination  of  fastener  loads  and  skin 
stress  gradients  can  be  accomplished  using  the  FRANC2D/L  program  developed 
jointly  at  Cornell  University  and  Kansas  State  University  sponsored  by  NASA- 
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Langley  Research  center.  The  finite  element  model  can  be  created  using  the- 
FRANC2D/L  preprocessor  program  called  CASCA  developed  at  Cornell 
University.  Figure  A2  shows  a  typical  finite  element  model  representing  one 
quarter  of  the  modified  skin  and  a  mounting  plate  attached  to  the  skin  created 
using  the  CASCA  program.  In  the  figure,  the  triangle  symbol  indicates  the 
fastener  location  and  the  heavy  solid  line  represents  the  boundary  of  the 
mounting  plate.  In  the  analysis,  the  skin  is  subjected  to  a  uniform  far  field  stress 
applied  along  the  top  boundary.  Due  to  double  symmetry  of  the  geometry  as  well 
as  the  loading,  translational  displacements  are  constrained  along  the  left  and 
bottom  edges  of  the  model  (both  skin  and  mounting  plate)  in  the  analysis. 


Figure  A2.  Finite  Element  of  One  Quarter  of  the  Modified  Skin 
with  Mounting  Plate  and  Fasteners 

In  FRANC2D/L  analysis,  the  fastener  shear  rigidity  needs  to  be  calculated 
separately  and  input  into  the  program.  The  FRANC2D/L  assembles  the  stiffness 
contributing  from  each  of  the  structural  components,  skin,  mounting  plate  and 
fasteners,  to  form  an  overall  stiffness  matrix.  With  the  loads  applied  to  the 
structural  model,  FRANC2D/L  solves  for  deformations,  load  transfers,  and  stress 
in  the  skin  and  mounting  plate. 

The  shear  rigidity  of  the  fastener  can  be  calculated  using  the  Swift’s  empirical 
equation: 

K  = _ Ed _ 

A  +  B(— + - - - ) 

^Skin  ^MountingPlate 
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where  E  is  the  weighted  Young’s  modules  of  the  skin  and  mounting  plate,  d  is  the- 
fastener  hole  diameter,  tswn  and  tMounting  Plate  are  the  thickness  of  the  skin  and 
mounting  plate,  respectively,  and  A  and  B  are  material  dependent  empirical 
constants.  The  values  of  A  and  B,  for  aluminum  fasteners,  are  5.0  and  0.8, 
respectively. 

In  the  analysis,  a  reference  uniform  far  field  stress  of  1,000  psi  is  applied  along 
the  top  edge  of  the  model. 

A4.  Parameters 

Parameters  considered  in  the  fastener  loads  and  skin  stress  gradients 
calculations  for  square/rectangular  and  circular  installations  include: 

Skin  thickness:  0.032",  0.040”,  0.050”,  and  0.063” 

Mounting  plate  thickness:  1  gauge  higher  than  skin  thickness 
Antenna  connector  hole  diameter:  1 .0”,  2.0”,  and  3.0” 

Fastener  size:  1/8”,  5/32”,  and  3/16” 

Number  of  fastener  rows:  1  ~  5  on  each  side 
Aspect  ratio:  1 :1  (square),  1 :1 .3  and  1 :1 .6 
Materials:  Aluminum  skin  and  mounting  plate 
Aluminum  fastener 

In  the  database  development,  the  number  of  fastener  rows,  in  rectangular 
installations,  normal  to  the  stress  application  direction  is  limited  to  3  on  each 
side.  However,  it  is  limited  to  5  on  each  side  for  the  fastener  rows  parallel  to  the 
stress  application  direction.  For  circular  installations,  the  number  of  fastener 
rings  is  limited  to  3. 

A  total  of  180  rectangular  and  108  circular  antenna  installations  are  included  in 
the  fastener  loads  and  skin  stress  gradients  database. 

A5.  Results 

Analysis  results  of  fastener  loads  along  the  critical  outmost  faster  row  are 
obtained  based  on  a  1,000  psi  reference  skin  stress  for  180  rectangular  and  108 
circular  antenna  installations.  These  fastener  loads  are  used  to  calculate  the 
bearing  and  bypass  stresses  which  are  used  in  the  stress  intensity  factor 
calculation.  Stress  gradients  along  a  crack  growth  path  normal  to  the  stress 
application  direction  are  also  obtained  for  a  crack  initiating  from  the  antenna 
connector  hole.  These  stress  gradients  are  used  in  the  stress  intensity  factor 
calculation. 
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A6.  Examples 


Figure  A3  depicts  an  example  problem  of  an  antenna  installation  on  the  fuselage 
of  a  commuter  airplane. 


Skin  thickness:  0.050” 

Mounting  plate  thickness:  0.063" 
Cutout  diameter:  1 .0" 

Fastener  diameter:  3/1 6” 

Fastener  rows:  3  on  each  side 
Fastener  pitch:  1.0” 

Skin  and  mounting  plate:  AI2024T3 
Fastener:  NASI  097 


Figure  A3.  Description  of  the  Example  Problem 


Fastener  loads  in  pounds  obtained  from  FRANC2DL  are  shown  in  Figure  A4.  It 
shows  that  load  transfers  at  fasteners  along  the  outermost  row  are  larger  than 
other  rows. 
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Figure  A4.  Fastener  Loads 
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One  of  the  initial  flaw  assumptions  in  an  antenna  installation  is  to  postulate  an 
initial  crack  at  the  antenna  connector  hole.  With  this  assumption,  the  crack 
continues  to  grow  along  a  path  normal  to  the  stress  application  direction  as 
shown  in  Figure  A5. 
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Figure  A5.  Potential  Crack  Growth  Paths 

The  stress  gradient  along  the  potential  crack  growth  path  is  shown  in  Figure  A6. 
The  distance  d  shown  in  the  figure  is  measured  from  the  center  of  the  hole.  This 
stress  gradient  is  used  in  the  stress  intensity  factor  calculation  for  damage 
tolerance  analysis. 


Figure  A6.  Stress  Gradient  along  Potential  Crack  Growth  Path 
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Appendix  B 

Load/Stress  Spectrum  Development 


B1.  Introduction 

This  Appendix  describes  procedures  used  in  RAPIDC  to  develop  the  load  spectrum  at 
the  airplane  C.G.  location  for  commuter  airplanes.  The  load  spectrum  consists  of  load 
sequences  that  an  airplane  experiences  in  a  flight  including  pre-flight  taxi,  climbing, 
cruise,  descending,  landing  impact,  and  post-flight  taxi.  To  perform  the  damage 
tolerance  analysis  of  a  structural  component  on  an  airplane,  a  stress  spectrum  is 
required.  In  RAPIDC,  a  simple  method  based  on  the  material  allowable  is  used  to 
derive  the  stress  spectrum  from  the  load  spectrum. 

The  steps  taken  to  generate  load  and  stress  spectra  in  RAPIDC  are  outlined  below. 

1  User  provides  RAPIDC  with  characteristics  of  the  airplane  such  as  the  design 
maximum  take-off  weight,  the  wing  surface  area,  the  maximum  design  cruise 
altitude,  and  the  normal  maximum  working  pressure. 

2  User  provides  RAPIDC  with  up  to  six  airplane  usage  data  such  as  the  flight 
distances,  the  cruise  altitude,  and  the  cruise  speed. 

3  Incorporate  maneuver  and  gust  exceedance  data  from  DOT/FAA/CT-91/20 
(Reference  1)  for  single  and  twin  engine  pressurized  general  usage  into  RAPIDC. 

4  Randomize  flight  profile  mixes  to  create  load  sequences  for  each  flight  segment  and 
form  a  complete  flight  load  spectrum. 

5  Derive  the  stress  spectra  at  the  location  of  interest  from  the  load  spectrum  based  on 
the  skin  material  allowable. 

The  load  spectrum  created  for  RAPIDC  consists  of  3,000  flights  of  load  sequences  in 
terms  of  the  change  of  airplane  C.G.  normal  acceleration  (Ag)  from  the  1-G  flight. 
RAPIDC  creates  a  load  sequence  of  one  tenth  of  the  airplane  design  life  when  the  user 
optionally  provides  RAPIDC  with  the  design  life  of  the  airplane. 
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B2.  Input  Parameters  Required  for  Load  and  Stress  Spectra  development 

To  generate  the  load  spectrum,  the  user  is  required  to  provide  RAPIDC  with  the 
following  data: 

1 .  Airplane  design  maximum  take-off  weight  (lbs) 

2.  Nominal  maximum  working  pressure  (psi) 

3.  Flight  distance  (statue  or  nautical  miles) 

4.  Maximum  design  cruise  altitude  (feet) 

5.  Cruise  speed  (KEAS) 

6.  Cruise  altitude  (feet) 

7.  Wing  surface  area  (feet2) 

8.  Airplane  design  life  (flights) 

The  following  additional  pieces  of  information  are  also  required  to  derive  the  stress 
spectrum  from  the  load  spectrum. 

1 .  Fuselage  radius  at  the  location  of  interest 

2.  Zone  number  indicating  the  location  of  interest  (see  Figure  B1) 

3.  Zone  Length  (S) 

4.  Distance  between  zone  starting  point  and  the  location  of  interest  (L) 

5.  Distance  between  cabin  floor  and  the  location  of  interest  (Z) 


Floor  Line 


Z  (negative  below  floor  line) 


Figure  B1 .  Airplane  Zones 

B3.  Flight  Profiles 

A  typical  flight  profile  of  commuter  airplanes  can  be  divided  into  stages  of  operation. 
The  time  line  of  the  operation  stages  is  schematically  shown  in  Figure  B2  and  the 
stages  used  in  the  RAPIDC  program  and  the  load  environments  are  described  in  Table 
B1.  RAPIDC  accepts  up  to  six  flight  profiles  provided  by  the  user.  A  randomized 
mission  mix  scheme  is  used  to  generate  the  load  spectrum. 
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Figure  B2.  Typical  Flight  Profile 


Table  B1 .  Flight  Segments 


Stage  No 


Stage 


Pre-flight  Taxi 
Climb 
Cruise 
Descent 
Post-flight  Taxi 


Load  Environment 


Ground 

Gust  and  Maneuver 
Gust  and  Maneuver 
Gust  and  Maneuver 
Ground 


B3.1  Assumptions 

To  construct  a  flight  profile,  the  following  assumptions  are  made: 

1  There  are  six  climbing  stages  and  six  descending  stages  for  a  flight  profile. 

2  For  a  flight  profile,  the  duration  for  climb  and  descent  is  each  determined  linearly 
using  three  baseline  flight  duration  data. 


3  The  airplane  equivalent  speed  increases  linearly  in  proportion  to  the  altitude  from 
the  time  of  takeoff  (zero  speed)  to  the  beginning  of  cruise.  Similarly,  it 
decreases  linearly  from  cruise  to  landing  (zero  speed)  during  descent. 
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The  three  baseline  flight  duration  data  in  RAPIDC  are  defined  as  follows. 


Baseline  1 : 


Flight  duration  =  30  minutes 

Climb  and  descent  duration  =  6  minutes  each 

Baseline  2: 


Flight  duration  =  60  minutes 

Climb  and  descent  duration  =  8  minutes  each 

Baseline  3: 


Flight  duration  =  120  minutes 

Climb  and  descent  duration  =  10  minutes  each 

B3.2  Construction  of  Flight  Profile 

The  following  steps  are  used  to  construct  a  flight  profile: 

1  Compute  the  altitude  for  the  climbing  stages  which  is  linearly  interpolated  based 
on  the  altitude  at  the  ground  zero  level  and  that  at  the  cruise  altitude. 

2  Compute  the  altitude  for  the  descending  stages  which  is  linearly  interpolated 
based  on  the  altitude  at  the  cruise  altitude  and  ground  zero  level. 

3  Compute  average  speed  for  the  climbing  stages  which  is  linearly  interpolated 
based  on  the  equivalent  speed  at  takeoff  and  that  at  cruise  as  a  function  of 
airplane  altitude. 

4  Compute  average  speed  for  the  descending  stages  which  is  linearly  interpolated 
based  on  the  equivalent  speed  at  cruise  and  that  at  landing  as  a  function  of 
airplane  altitude. 

5  Compute  the  duration  for  the  climbing  stages  which  is  linearly  interpolated  using 
three  baseline  flight  duration  data. 

6  Compute  the  duration  for  the  descending  stages  which  is  linearly  interpolated 
using  three  baseline  flight  duration  data. 

7  Compute  the  cabin  pressure  at  each  stage  as  a  function  of  altitude  from  the 
pressure  schedule  described  in  Section  B4.1.4. 
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B3.3  An  Example:  Flight  Profile 

Table  2  gives  parameters  used  to  create  a  flight  profile  for  a  typical  commuter  airplane 


Table  B2.  Operation  Parameter  Profile  for  a  Typical  Commuter  Airplane 


Airplane  Parameters 

Units 

Typical  Commuter 
Airplane 

Flight  Distance 

st.  miles 

204 

Design  Maximum 
Takeoff  Weight 

lb. 

16300 

Wing  Area 

ft2 

310 

Cruise  Altitude 

ft. 

35,000 

Maximum  Cruise 
Altitude 

ft. 

45000 

Nominal  Maximum 
Working  Pressure 

psi 

9.3 

Cruise  Speed 

KEAS 

240 

Climb  Duration 

hrs. 

0.12834 

Descent  Duration 

hrs. 

0.12834 

The  flight  profile  generated  from  above  operation  parameters  is  shown  in  Table  B3. 
Table  3.  A  Typical  flight  Profile  for  Commuter  airplanes 


Segment 

Altitude 

Speed 

Duration 

Number 

(1000  Ft.) 

(KEAS) 

(Hours) 

(psi) 

1 

^  0 

0.0 

0.0 

0.60 

2 

0-5.83 

20.0 

0.02139 

2.04 

3 

5.83-11.67 

60.0 

0.02139 

4.43 

4 

11.67-17.50 

100.0 

0.02139 

6.28 

5 

17.50-23.33 

140.0 

0.02139 

7.67 

6 

23.33-29.17 

180.0 

0.02139" 

8.68 

7 

29.17-35.0 

220.0 

0.02139 

9.35 

8 

35 

240.0 

0.59330 

9.59 

9 

35-29.17 

220.0 

0.02139 

9.35 

10 

29.17-23.33 

180.0 

0.02139 

8.68 

11 

23.33-17.50 

140.0 

0.02139 

7.67 

12 

17.50-11.67 

100.0 

0.02139 

6.28 

13 

1 1 .67-5.83 

60.0 

0.02139 

4.43 

14 

5.83-0 

20.0 

0.02139 

2.04 

15 

0 

0.0 

0.0 

0.60 

B4.  Load  Spectrum  at  Airplane  Center-of-Gravity 

This  section  describes  procedures  of  load  sequence  generation  at  the  airplane  C.G. 
location. 

B4.1  Load  Environments 

Load  environments  contributing  to  the  crack  growth  of  a  modified  skin  include  the 
following: 

1  In-flight  maneuver  operations  during  climb,  cruise,  descent  and  landing  impact 

2  Gust  loads  due  to  air  turbulence 

3  Ground  operations  such  as  pre-flight  taxi  and  post-flight  taxi 

4  Cabin  pressure 

The  frequency  distribution  of  the  accelerations  that  the  airplane  experiences  at  the 
center  of  gravity  can  be  expressed  in  tabulated  forms  such  as  the  load  exceedance 
tables.  These  exceedance  data  are  obtained  from  the  FAA  document  (Reference  1)  for 
general  aviation  aircraft.  The  acceleration  exceedance  data  for  each  load  environment 
are  presented  as  follows. 

B4.1.1  In-flight  Maneuver 

The  distributions  of  airplane  acceleration  during  climb,  cruise,  descent  and  landing 
impact  for  single  and  twin-engine  pressurized  general  usage  airplanes  are  shown  in 
Figure  B3.  Shown  in  the  figure  are  the  cumulative  frequency  of  exceedance  or 
cumulative  number  of  occurrences  per.  nautical  mile  as  a  function  of  acceleration 
fraction  an/an,LLF  where,  an  is  the  change  of  normal  acceleration  in  the  unit  of  g’s  and 
an,LLFis  the  change  of  limit  load  factor. 


■1  -0.75  -0.5  -0.25  0  0.25  0.5  0.75 

Acceleration  Fraction  a*  /  a*.uj = 


Figure  B3.  Maneuver  Exceedance  Data 
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It  is  noted  that  the  above  exceedance  data  are  the  cumulative  number  of  occurrences 
per  nautical  mile.  Load  spectra,  however,  are  generally  expressed  in  terms  of  cycles 
per  1,000  flights.  For  this  reason,  the  exceedance  or  the  cumulative  number  of 
occurrences  per  nautical  mile  needs  to  be  converted  to  that  per  1 ,000  flights.  This 
conversion  is  accomplished  by  multiplying  the  exceedance  per  nautical  mile  by  the  flight 
distance  in  nautical  miles  times  1,000  for  1,000  flights.  To  obtain  the  exceedance  per 
spectrum,  multiply  the  exceedance  per  1,000  flights  by  the  number  of  flights  per 
spectrum  (assuming  1/10  of  the  airplane  design  life,  default  to  3,000  flights)  and  then 
divided  by  1 ,000  flights. 

It  should  also  be  noted  that  the  exceedance  obtained  from  Reference  1  is  tabulated  in 
terms  of  acceleration  fraction  a n/  an,i_LF.  For  practical  use,  the  exceedance  is  expressed 
as  a  function  of  the  change  of  normal  acceleration  an  from  the  1  -G  stress  level.  The 
conversion  from  the  acceleration  fraction  a n/  an,LLF  to  the  change  of  normal  acceleration 
a„  can  be  accomplished  as  follows. 

For  positive  maneuver,  the  limit  load  factor  ni  is  calculated  as 

?4-000- 
W  +  10,000 

where  W  is  the  design  maximum  take-off  weight  of  the  airplane  in  pounds.  The 
calculated  value  of  ni  should  be  bounded  between  2.5  and  3.8.  The  change  of  limit 
load  factor  an,LLF  is  then  calculated  as 

an,UF  ~  ni  —  TO 

For  negative  maneuver,  the  limit  load  factor  nz  is  calculated  as 

n2  =  -0.5  n{ 

The  change  of  limit  load  factor  an,LiF  is  then  calculated  as 

an,LLF  =n2-\.0 

The  change  of  acceleration  an  is  then  obtained  by  multiplying  the  acceleration  fraction 
an  /  an,LLF  by  the  change  of  limit  load  factor  an,i_LF. 

Figure  B4  shows  an  example  of  the  maneuver  exceedance  per  1 ,000  flights  as  a 
function  of  the  change  of  normal  acceleration  for  the  operation  parameters  described  in 
Table  B2. 


B-7 


1.00E+07 
1  1.00E+06 
j=  1.00E+05 
o  1.00E+04 
8  1.00E+03 
^  1.00E+02 
o  1.00E+01 
■g  1.00E+00 
g  1.00E-01 
ft  1.00E-02 
1.00E-03 

-2.5  -2  -1.5  -1  -0.5  0  0.5  1  1.5  2  2.5 

Change  of  Normal  Acceleration  an 

Figure  B4.  An  Example  of  Maneuver  Exceedance 

B4.1.4  Gust 

The  distributions  of  airplane  acceleration  due  to  gust  for  single  and  twin-engine 
pressurized  general  usage  airplanes  are  shown  in  Figure  B5.  Shown  in  the  figure  are 
the  cumulative  number  of  occurrences  per  nautical  mile  as  a  function  of  acceleration 
fraction  an/an,LLF  where,  an  is  the  change  of  normal  acceleration  in  the  unit  of  g’s  and 
3n,LLF  is  the  change  of  limit  load  factor. 


Acceleration  Fraction  a*  /  a^up 

Figure  B5.  Gust  Exceedance  Data 
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The  gust  exceedance  data  obtained  from  Reference  1  are  the  cumulative  occurrences 
per  nautical  mile  in  terms  of  the  acceleration  fraction.  For  the  same  reasons  addressed 
in  Section  B4.1.2,  the  exceedance  needs  to  be  converted  to  cumulative  occurrences 
per  1 ,000  flights  as  a  function  of  the  change  of  normal  acceleration  from  the  1  -G  stress 
level.  This  can  be  accomplished  using  the  same  approach  described  in  Section  B4.1 .2. 
However,  the  positive  and  negative  gust  limit  load  factors  n3  and  n4  as  shown  in  Figures 
B6  and  B7,  respectively,  should  be  used  per  FAR  Part  23,  Appendix  A  (Reference  2). 


Figure  B6.  Determination  of  03  Factor  at  Speed  Vc 


"is 


Figure  B7.  Determination  of  n4  Factor  at  Speed  Vc 

In  Figures  B6  and  B7,  nt  is  the  positive  maneuver  limit  load  factor,  W  is  the  design 
maximum  take-off  weight  of  the  airplane  (pounds),  S  is  the  airplane  wing  surface  (feet ), 
Vc,  sei  is  the  airplane  cruise  speed  (KTAS),  and  Vc,Min  is  the  minimum  design  cruise 
speed  (KTAS)  which  is  calculated  as 


W 

^CMin  ~  T7.0y«t  — 
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Once  the  positive  and  negative  limit  load  factors  n3and  n4are  obtained,  the  change  of 
limit  load  factors  can  be  calculated  as 

For  positive  gust  exceedance 


an,LLF  =«3-L0 

For  negative  gust  exceedance 

an,LLF  =«4-1*0 

Figure  B8  shows  an  example  of  the  gust  exceedance  per  1,000  flights  as  a  function  of 
the  change  of  normal  acceleration  for  the  operation  parameters  described  in  Table  B2. 


1.00E+07 
1.00E+06 
|  1.00E+05 
j=  1.00E+04 
§  1.00E+03 
1.00E+02 
o  1.00E+01 
a  1.00E+00 
|  1.00E-01 
i£  1.00E-02 
1.00E-03 
1.00E-04 

-3-2-10  1  2  3 

Change  of  Normal  Acceleration  an 

Figure  B8.  An  Example  of  Gust  Exceedance 


B4.1.3  Taxi 

The  distributions  of  airplane  acceleration  during  ground  operations  are  shown  in  Figure  - 
B9  for  1 ,000  flight  operations. 
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Figure  B9.  Exceedance  Data  of  Taxi 


B4.1.4  Differential  Cabin  Pressure 

The  pressure  differential  AP  at  an  airplane  altitude  Hpiane  is  the  pressure  difference 
between  the  cabin  pressure  Pcawn  and  the  atmospheric  pressure  Patm.  The  cabin 
pressure  is  generally  expressed  in  terms  of  the  cabin  altitude  Hcabin,  which  is  limited  to 
8,000  feet  under  normal  operating  conditions  as  set  forth  in  FAR  25.841 .  The 
atmospheric  pressure  at  an  altitude  can  be  obtained  from  the  Standard  Atmosphere 
table  published  by  ICAO  International  Civil  Aviation  Organization  (Reference  3).  It  is 
assumed  that  the  cabin  altitude  Hcabin  is  related  to  the  airplane  altitude  Hpiane  in  a 
parabolic  form  as  described  below: 


Hcabin=  Cp  (Hpi3ne) 


where  Hcabin  and  Hpiane  are  the  cabin  and  airplane  altitudes  in  feet,  respectively,  n  is  the 
cabin  pressurization  factor,  and  Cp  is  a  coefficient  calibrating  the  nominal  maximum 
working  pressure  at  the  maximum  airplane  altitude.  The  cabin  pressurization  factor  n 
takes  a  value  of  1 .6  in  RAPIDC. 

To  determine  the  pressure  differential  Ap  in  psi  at  an  airplane  altitude  in  feet,  the  user 
provides  RAPIDC  with,  the  maximum  airplane  altitude  Hmax  in  feet  and  the  nominal 
maximum  working  pressure  Pnominai.max  in  psi.  Using  the  steps  as  described  below, 
RAPIDC  calculates  the  pressure  differential  as  a  function  of  cruise  atitude.  In  the 
calculation,  the  altitude  and  the  pressure  are  in  the  units  of  feet  and  psi,  respectively. 
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(1)  Calculate  the  cabin  altitude  Hcabin 

(a)  Determine  the  maximum  pressure  pH, max  at  the  maximum  airplane 
altitude  Hmax 

(b)  Add  the  nominal  maximum  pressure  Pnominai.max  to  pH.max  to  obtain  the 
cabin  pressure  pcabin 

(c)  Determine  the  cabin  altitude  in  feet  Hcabin  at  cabin  pressure  pCabin 

(2)  Calculate  the  coefficient  Cp  using  the  equation 

M 

r*  —  11  cabin 


(3)  Calculate  the  pressure  differential  Ap 

(a)  Determine  the  atmospheric  pressure  Patm  at  the  airplane  altitude  Hpiane 

(b)  Calculate  the  cabin  altitude  Hcabin  and  determine  the  cabin  pressure 

Pcabin 

(c)  Calculate  the  pressure  differential  Ap  =  pCabin  -  Patm 

(d)  Add  0.6  psi  to  Ap  per  FAR  23.574 

As  an  example,  shown  in  Figure  BIO  below  is  the  pressure  schedule  for  an  airplane 
providing  a  9.3  psi  nominal  maximum  working  pressure  at  45,000  feet  airplane 
maximum  altitude. 


Figure  BIO.  Example  of  Pressure  Schedule 

B4.2  Creation  of  Load  Sequences 

The  load  spectra  are  created  according  to  the  typical  airplane  operation  profile 
described  in  Table  B1  for  commuters.  The  procedures  of  load  spectrum  generation 
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are  coded  in  the  FORTRAN  computer  program  "CGENLINT"  and  are  briefly  described 
below: 

1  Create  load  exceedance  tables 

The  exceedance  tables  are  calculated  at  0.1  g  increments  starting  at  O.Og  for  a 
user  specified  number  of  flights  for  the  gust  and  maneuver  load  environments. 
The  exceedances  are  converted  from  nautical  miles  to  exceedances  per  1 ,000 
flights  based  on  the  flight  distance. 

2  Create  load  occurrence  tables 

The  occurrence  tables  are  calculated  at  0.1  g  increments  starting  at  0.1 5g  based 
on  the  exceedances.  The  gust  and  maneuver  occurrence  tables  per  1,000 
flights  are  converted  to  occurrence  tables  per  spectrum  based  on  the  user 
specified  number  of  flights.  Similarly,  the  taxi  occurrence  table  per  1,000 
landings  is  converted  to  occurrences  per  spectrum  based  on  the  user  specified 
number  of  flights.  For  each  load  environment  table,  the  total  number  of  positive 
occurrences  per  spectrum  and  the  total  number  of  negative  occurrences  per 
spectrum  should  be  equal.  If  any  difference  between  the  total  numbers,  it  is 
correspondingly  balanced  at  0.0  Ag.  The  occurrence  tables  are  also  referred  as 
“pool  of  load  cycles”. 

3  Create  load  sequence  for  each  flight  profile  stage 

The  load  sequence  is  generated  by  randomly  picking  cycles  from  the  appropriate 
pools  for  each  flight  profile  stage.  A  Ag  value  is  randomly  picked  from  the 
positive  pool  of  load  cycles,  and  paired  with  another  Ag  value  that  is  randomly 
picked  from  the  negative  pool  of  load  cycles.  Each  cycle  can  only  be  picked 
once  and  every  cycle  in  the  pools  is  used.  For  each  stage,  the  number  of  cycles 
to  be  picked  is  determined  by  the  duration  of  that  stage  and  the  number  of  cycles 
in  the  pool  as  follows: 


Npick=  integer  { Npooi  x  (Tstage/Tpooi) } +  (1)nst 

where 

Npick  =  number  of  cycles  to  be  picked  for  the  stage 

Npooi  =  number  of  cycles  in  the  original  pool 

Tstage  =duration  of  the  stage  for  one  flight 

Tpooi  =  duration  of  the  pool  for  a  user  specified  number  of  flights 
+  (1)nSf=  additional  cycle  that  is  added  to  the  flight  every  time  the  accumulation  of 
the  fraction  portion  of  {Np00i  x  (Tstag</Tpooi)}  becomes  one  cycle.  The 
number  of  flights  that  have  one  additional  cycle  is  determined  as  follows: 

nsf=  NPooi  -  integer{Np00i  x  (Tstage/Tpooi)}  x  a  user  specified  number  of  flights 
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In  general,  the  cyclic  loads  are  paired  with  respect  to  the  1-g  state  to  form  a  full  cycle. 
The  order  of  load  sequences  for  airborne  stages  are  as  follows:  gust  and  maneuver. 
The  minimum  and  maximum  load  factors  and  number  of  cycles  for  each  operation 
environment  are  discussed  in  section  B5. 

B5.  Contents  of  Load  Spectrum 

The  generic  load  spectra  are  stored  in  binary  form  and  can  be  retrieved  using  the 
following  FORTRAN  code: 

READ {iu)  AN1  ,NSEG,(CPD(I),I=1  ,NSEG) 

READ(/u)  NFLTS 
DO  1=1, NFLTS 

READ(iu)  IFLT,IPAIR,((IDELTG(K,J),K=1 ,2),J=1  ,IPAIR), 

+  ((ICYC(K,J),K=1  ,NSEG),J=1 ,2) 

END  DO 


where 

iu  is  the  FORTRAN  unit  of  the  spectrum  file 

AN1  is  the  positive  maneuver  limit  load  factor 

NSEG  is  the  number  of  stages  in  the  flight  profile 

CPD  is  the  cabin  pressure  differentials  for  each  stage  of  the  flight  profile; 
REAL*8 

NFLTS  is  the  total  number  of  flights  in  the  spectrum  (i.e.,  default  to  3,000) 

IFLT  is  the  flight  number  (i.e.,  1  through  3,000) 

IPAIR  is  the  number  of  incremental  load  factors  in  the  flight  (max=2,000) 

IDELTG  is  the  array  of  incremental  load  factors;  to  save  on  storage,  the 
incremental  load  factors  are  multiplied  by  1000,  truncated,  and  stored  in 
the  INTEGER*2  array 

ICYC  is  the  array  of  the  summation  count  of  the  number  of  incremental  load 
factors  for  each  stage  of  the  flight  profile 

B6.  Derivation  of  Stress  Spectrum 

The  stress  spectrum  at  the  location  where  the  damage  tolerance  analysis  is  performed, 
in  either  the  longitudinal  or  circumferential  direction,  can  be  derived  from  the  load 
spectrum  using  the  following  equation: 

a  =  ap  +  (l  ±  Ag)  dig 
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Where  cr  is  the  min.  or  max.  stress  of  a  stress  cycle  in  the  stress  spectrum 
Op  is  the  stress  due  to  cabin  pressure  at  the  location 

Gig  is  the  1-g  stress  due  to  airplane  inertia  and  aerodynamic  loads  at  the  location 

A g  is  the  change  of  normal  acceleration  in  the  unit  of  g’s  in  the  load  spectrum 

To  estimate  the  longitudinal  1-g  stress,  Gig  at  the  location  where  the  damage  tolerance 
analysis  is  performed,  the  following  assumptions  are  made. 

(1)  The  longitudinal  stress  in  flight  at  the  crown  location  can  be  described  as 
shown  in  Figure  B1 1.  The  stress  in  Zone  I  can  be  linearly  interpolated  within 
the  zone. 


pr/2t+  a  lgm(a 


pr/2t 

II 

III 

^  Lone  i  V 

Figure  B1 1  Longitudinal  Stress  in  Flight  at  the  Crown  Location 

(2)  The  longitudinal  stress  on  the  ground  at  the  crown  location  can  be  described 
as  shown  in  Figure  B12.  The  stress  in  Zone  II  can  be  linearly  interpolated 
within  the  zone. 
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Figure  B1 2.  Longitudinal  Stress  on  the  Ground  at  the  Crown  Location 

(3)  The  longitudinal  stress  between  the  crown  and  belly  locations  can  be  linearly 
interpolated  as  shown  in  Figure  B13. 


pr'2t+alt 


Figure  B1 3.  Longitudinal  Stress  at  a  Distance  Z  from  cabin  Floor 
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With  these  assumptions,  the  longitudinal  stress  at  a  location  can  be  estimated  once  the 
maximum  design  1-G  stress,  crig,  max  is  determined.  The  determination  of  the  cr1g,  max  is 
described  as  follows. 


B6.1  Determination  of  Maximum  Design  1-G  Stress 


It  is  assumed  that  the  skin  stress  due  to  pressure  and  the  maximum  design  1-g  stress, 

c 3ig,max  at  the  limit  load  is  equal  to  the  design  limit  stress  of  the  fuselage  skin.  This 
assumption  can  be  expressed  in  an  equation  as 


O  +  0.6)£ 
It 


+  «,  G 


lg’.max 


where  p  is  the  pressure  differential,  ni  is  the  maneuvering  limit  load  factor,  Ftu  is  the 
ultimate  tensile  strength  of  the  skin  material,  C  is  the  material  allowable  reduction  factor 
for  assembled  structure,  and  R  and  t  are  fuselage  radius  and  skin  thickness  at  the 
location,  respectively.  The  additional  pressure  0.6  psi  added  to  the  pressure  differential 
is  to  account  for  the  external  aerodynamic  suction  pressure  (0.5  psi)  plus  regulator 
tolerance  (0.1  psi)  per  FAR  23.574.  In  the  above  equation,  n1  is  the  positive  maneuver 
limit  load  factor  which  can  be  calculated  using  the  equation  described  in  Section  B4.1.1. 
An  empirical  value  of  0.88  is  used  for  the  material  allowable  reduction  factor,  i.e., 
C=0.88. 


Using  the  above  equation,  the  maximum  design  1-g  stress,  Gig,max  at  the  limit  load  can 
be  determined  as 


<7 


l£,max 


(p  +  0.6)  /?/ 


"i 


For  ground  operation  conditions,  the  maximum  1-g  stress  is  assumed  to  be  1/3  of 
Gig, max  as  shown  in  Figure  B12. 

B6.2  One-G  Stress  at  Location  of  Interest 


Once  the  maximum  design  1-g  stress,  Gig,max  at  the  limit  load  is  determined,  the  1-g 
stress,  <7ig  at  a  location  of  interest  can  be  determined  as  follows. 

(A)  In-flight  conditions 


(i)  Longitudinal  stress  in  Zone 


<yig=\ 


(-) 

UJ 

UJ 

l£,max 
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(ii)  Longitudinal  stress  in  Zones  II,  III,  and  IV: 


kRj 


1  !g,max 


iii)  Circumferential  stress  in  all  Zones: 


°rtg  =0 


(B)  Ground  conditions 

(i)  Longitudinal  stress  in  Zone  I: 


g1s=0 


(ii)  Longitudinal  stress  in  Zones  II: 


.  1 

ai*s*7 


(L) 

UJ 

UJ 

l£,max 


iii)  Longitudinal  stress  in  Zones  III  and  IV 


.  1 
=  A,  - 


(Z) 

—  R 

UJ 


(iv)  Circumferential  stress  in  all  Zones 


als=0 


where 


Z  =  the  distance  between  the  location  and  the  cabin  floor  line,  positive  for 
locations  above  the  floor  line;  negative  for  locations  below  the  floor  line 

R  =  the  radius  of  the  fuselage 

L  =  the  distance  between  the  location  and  the  reference  point  in  the  zone  (see 
Figure  B1) 

S  =  the  length  of  zone  (see  Figure  B1 ) 

A,  =  assumed  average  payload  factor  equal  to  0.70 
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B7.  Rain-flow  Count  and  Truncation 


The  stress  sequences  at  the  location  is  rain-fiow  counted  for  each  flight  and  is 
truncated  at  the  range  of  2000  psi.  The  rain-flow  counting  method  is  one  of  the  re¬ 
sequence  processes  generally  applied  to  a  flight-by-flight  stress  spectrum  prior  to  the 
crack  growth  analysis.  The  method  is  illustrated  in  figure  B13.  Orienting  the  graphical 
display  of  the  stress  sequence  vertically,  the  stress  spectrum  is  regarded  as  a  stack  of 
tiles  on  the  roof.  Rain  continues  to  flow  from  one  tile  to  the  next.  If  the  water  runs  off 
the  tile,  it  drops  down  on  the  tile  below.  It  is  assumed  in  this  process  that  when  the  tile 
is  already  wet,  the  water  stops  to  flow.  Using  this  analogy,  the  stress  range  is  regarded 
as  the  rain-flow  range  indicated  by  AB,  CD,  etc.,  in  Figure  B13.  The  rain-flow  counting 
process  has  been  standardized  and  is  described  in  ASTM  E-1049  (Reference  4). 


Figure  B13.  Rain-flow  Count 

The  truncation  of  lower  stress  levels  in  a  stress  spectrum  is  employed  in  RAPIDC  for 
efficiency  in  crack  growth  analysis.  The  stress  truncation  omits  stress  cycles  below  a 
certain  magnitude.  It  is  a  general  practice  to  accelerate  the  calculation.  The  argument 
is  that  low  stress  excursions  do  not  contribute  much  to  crack  growth,  especially  in  view 
of  the  retardation  effects.  The  final  stress  sequences,  after  rain-flow  counting  and 
truncation,  are  used  in  the  crack  growth  analysis. 
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Appendix  C 

One-Cycle  Equivalent  Stress  Calculation 


Cl.  Introduction 

This  appendix  describes  the  derivation  of  the  one-cycle  equivalent  stress  for  a  stress 
spectrum  either  generated  by  RAPIDC  or  provided  by  the  user.  The  one-cycle 
equivalent  stress  is  used  in  the  crack  growth  prediction  of  the  modified  skin  to  antenna 
installation  using  the  simplified  method. 

C2.  Derivation 

in  the  crack  growth  analysis  using  the  simplified  method,  the  crack  growth  life  Ng  for  a 
crack  growing  from  the  size  ^  to  the  size  a.  (a.  >  a.)  under  a  repeated  flight  can  be 

calculated  using  the  following  equation  [1]: 

Nir|(SGijrp  “  0) 

In  equation  1 ,  C  and  p  are  the  coefficients  in  the  Walker's  crack  growth  equation, 

^  =  C{(1-R)qKMax}P  (2) 

where  Kmax  is  the  stress-intensity  factor  evaluated  at  the  stress  aMax  of  the  given  stress 
cycle  and  R  is  the  stress  ratio  of  that  cycle. 

The  parameter  S  in  equation  1  is  the  one-cycle  equivalent  stress  of  the  repeated  flight, 
and  Gjj  is  the  geometry  term  pertinent  to  the  crack  geometry  of  the  modified  skin. 

The  equivalent  stress  S  can  be  obtained  by  the  equation: 

s={iKxJ<i-Ri>q]> p  <3> 

j=i 

where  q  is  the  coefficient  in  the  Walker’s  crack  growth  equation,  crMax  j  and  Rj  are  the  k 

maximum  stress  and  the  stress  ratio,  respectively,  of  the  j-th  cycle,  and  m  is  the  total 
number  of  cycles  in  the  repeated  flight. 

The  geometry  term  Gjj  can  be  calculated  by  the  equation: 

Gq  =  {J*‘[p(a)Vita]_pdar1/p  (4) 
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where  p(a)  is  the  geometry  factor,  and  a  is  the  crack  length. 

Suppose  a  stress  spectrum  at  the  antenna  location  contains  M  successive  flights.  Let 
mj  be  the  number  of  stress  cycles  in  the  j-th  flight.  A  representative  stress  spectrum  is 
depicted  below  in  Figure  1 . 


Figure  1 .  A  Representative  Stress  Spectrum 

With  the  stress  spectrum  in  each  flight  prescribed,  the  one-cycle  equivalent  stress  Sj  for 
each  j-th  flight  can  be  calculated  using  equation  3,  and  is  schematically  shown  below  in 

Figure  2. 


Figure  2.  One-Cycle  Equivalent  Stress 

Consider  a  stress  spectrum  containing  M  flights  as  shown  in  Figure  2.  The  equivalent 
stress  SEq  of  a  flight  in  the  spectrum,  as  shown  below  in  Figure  3 


Figure  3.  Equivalent  Stress  of  M  Flights 
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can  be  calculated  as 

M 

zv 

1=1 

M 


Substitution  of  equation  3  into  equation  5  gives 

M  f  mi 

IlW-vT 

i=1  Itl _ 

M 


which  can  be  rewritten  as 


The  variable  n  in  equation  7  is  the  total  number  of  stress  cycles  in  the  stress  spectrum. 

C3.  Implementation  and  Validation 

The  calculation  of  the  equivalent  stress  Seq  for  a  stress  spectrum  containing  M  flights 
has  been  implemented  in  RAPIDC.  The  implementation  was  previously  validated  in 
RAPID,  a  precedent  of  RAPIDC,  through  an  example.  In  order  not  to  repeat  the 
validation  process,  the  previous  example  using  a  skin  repair  is  reprinted  as  follows. 

A  skin  repair  with  an  external  doubler  mechanically  fastened  to  the  skin  is  shown  below 
in  Figure  4.  The  skin  and  repair  doublers  are  made  of  the  2024-T3  clad  sheet.  Two 
fastener  types,  NAS1097-E6  and  HL326-6  are  used.  Dimensions  of  the  skin  cutout  and 
the  repair  doublers  are  given  in  the  figure. 
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2024-T3  Clad  Sheet 
Thickness  *  0.063” 
Doubler 

2024-T3  Clad  Sheet 
Thickness  *  0.071" 
Rivets: 

0  NAS1097-E6 
•  HL326-6 
Pitch  -  1.0" 

Edge  distance  -  0.5" 


Figure  4.  Description  of  the  Example  Repair 

The  radius  near  the  repair  is  120.0  inches,  and  the  pressure  differential  is  assumed  to 
be  8.6  psi.  To  perform  the  crack  growth  analysis,  an  initial  longitudinal  skin  crack  of 
0.05”  is  postulated  at  the  critical  center  fastener  hole  in  the  first  fastener  row. 

It  is  assumed  that  the  skin  near  the  repair  location  is  subjected  to  a  stress  spectrum  of 
20,887  stress  cycles  representing  1/10  of  the  design  life  of  the  airplane.  An  overall 
spectrum  summary  is  provided  in  the  table  1  below. 
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Table  1  Overall  Spectrum  Summary 


PEAK 

DISTRIBUTION 

RANGE 

DISTRIBUTION 

STRESS  RATIO 
DISTRIBUTION 

Peak  Stress 

Cvcles 

Stress  Ranee 

Cvcles 

Stress  Ratio 

Cycles 

Below  -5000. 

0 

Below  0. 

Below  -2.00 

0 

-5000.  -4000. 

0 

0. 

1000. 

0 

-2.00  -1.90 

0 

-4000.  -3000. 

0 

1000. 

2000. 

0 

-1.90  -1.80 

0 

-3000  -2000. 

0 

2000. 

3000. 

14585 

-1.80  -1.70  ' 

0 

-2000.  -1000. 

0 

3000. 

4000. 

379 

-1.70  -1.60 

0 

-1000.  0. 

0 

4000. 

5000. 

13 

-1.60  -1.50 

0 

0.  1000. 

0 

5000. 

6000. 

0 

-1.50  -1.40 

0 

1000.  2000. 

0 

6000. 

7000. 

2361 

-1.40  -1.30 

0 

2000.  3000. 

0 

7000. 

8000. 

3256 

-1.30  -1.20 

0 

3000.  4000. 

0 

8000. 

9000. 

270 

-1.20  -1.10 

0 

4000.  5000. 

145 

9000. 

10000. 

22 

-1.10  -1.00 

0 

5000.  6000. 

7358 

10000. 

1000. 

1 

-1.00  -0.90 

0 

6000.  7000. 

1717 

11000. 

12000. 

0 

-0.90  -0.80 

0 

7000.  8000. 

802 

12000. 

13000. 

0 

-0.80  -0.70 

0 

8000.  9000. 

9697 

13000. 

14000. 

0 

-0.70  -0.60 

0 

9000.  10000. 

1085 

14000. 

15000. 

0 

-0.60  -0.50 

0 

10000.  11000. 

78 

15000. 

16000. 

0 

-0.50  -0.40 

0 

11000.  12000. 

5 

16000. 

17000. 

0 

-0.40  -0.30 

0 

12000.  13000. 

0 

17000. 

18000. 

0 

-0.30  -0.20 

0 

13000.  14000. 

0 

18000. 

19000. 

0 

-0.20  -0.10 

0 

14000.  15000. 

0 

19000. 

20000. 

0 

-0.10  0.00 

0 

15000.  16000. 

0 

20000. 

21000. 

0 

0.00  0.10 

0 

16000.  17000. 

0 

21000. 

22000. 

0 

0.10  0.20 

4571 

17000.  18000. 

0 

22000. 

23000. 

0 

0.20  0.30 

1340 

18000.  19000. 

0 

23000. 

24000. 

0 

0.30  0.40 

17 

19000.  20000. 

0 

24000. 

25000. 

0 

0.40  0.50 

126 

20000.  21000. 

0 

25000. 

26000. 

0 

0.50  0.60 

2122 

21000.  22000. 

0 

26000. 

27000. 

0 

0.60  0.70 

8453 

22000.  23000. 

0 

27000. 

28000. 

0 

0.70  0.80 

4258 

23000.  24000. 

0 

28000. 

29000. 

0 

0.80  0.90 

0 

24000.  25000. 

0 

29000. 

30000. 

0 

0.90  1.00 

0 

25000.  26000. 

0 

30000. 

31000. 

0 

1.00  1.10 

0 

26000.  27000. 

0 

31000. 

32000. 

0 

1.10  1.20 

0 

27000.  28000. 

0 

32000. 

33000. 

0 

1.20  1.30 

0 

28000.  29000. 

0 

33000. 

34000. 

0 

1.30  1.40 

0 

29000.  30000. 

0 

34000. 

35000. 

0 

1.40  1.50 

0 

30000.  31000. 

0 

35000. 

36000. 

0 

1.50  1.60 

0 

31000.  32000. 

0 

36000. 

37000. 

0 

1.60  1.70 

0 

32000.  33000. 

0 

37000. 

38000. 

0 

1.70  1.80 

0 

33000.  34000. 

0 

38000. 

39000. 

0 

1.80  1.90 

0 

34000.  35000. 

0 

39000. 

40000. 

0 

1.90  2.00 

0 

Above  35000. 

0 

Above  40000. 

0 

Above  0.0 

0 

Total  Cycles  =  20887 

Total  No.  of  Flights  =  5910 

Highest  Peak  =  11760.30 

Lowest  Valley  =  1082.95 

Largest  Range  =  10283.55 
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To  perform  the  crack  growth  analysis,  it  is  assumed  that  a  crack  of  size  equal  to  0:05 
initially  exists  at  the  critical  fastener  hole  in  the  skin.  The  crack  grows  towards  the 
adjacent  hole  and  enters  the  hole.  When  the  crack  grows  into  the  hole,  it  is  assumed 
that  two  equal-length  cracks  of  0.005"  exist  at  outer  holes.  The  crack  continues  to  grow 
to  enter  next  adjacent  holes.  The  analysis  proceeds  until  the  residual  strength  of  the 
repaired  skin  is  no  longer  greater  than  the  limit  stress  of  the  skin. 

The  scenarios  for  (a)  a  crack  growing  towards  the  adjacent  hole  and  (b)  continuing 
damage  of  the  crack  growing  into  the  hole  are  depicted  below  in  Figure  5. 


(a)  Prior  to  Entering  the  Hole 


O'Bypass 
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aQross 

(b)  Continuing  Damage 
Figure  5.  Crack  Growth  Scenario 

To  perform  the  crack  growth  analysis  using  the  simplified  method,  an  equivalent  stress 
Sea  equal  to  11.915  ksi  was  calculated  using  equation  7.  With  this  equivalent  stress, 
RAPID  calculates  the  crack  growth.  The  analysis  result  was  then  compared  with  that 
obtained  using  the  RAPID  cycle-by-cycle  method. 

Crack  growth  analysis  results  obtained  from  RAPID  using  the  equivalent  stress  in  the 
simplified  method  and  the  cycle-by-cycle  method,  are  plotted  in  Figure  6. 
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Figure  6.  Crack  Growth  Life  of  the  Example  Repair 

As  shown  in  Figure  6,  the  crack  growth  obtained  using  the  equivalent  stress  in  the 
simplified  method  is  nearly  identical  to  that  from  the  cycle-by-cycle  method.  In  this 
analysis,  the  crack  growth  was  calculated  based  on  the  zero  stress  ratio  (R  =  0.0)  in  the 
simplified  method.  The  crack  growth  da/dN  material  data  used  in  the  cycle-by-cycle 
method  were  obtained  for  various  stress  R  ratios  using  the  Walker’s  equation.  In  this 
manner,  the  straight-line  portion  of  the  da/dN  data  were  used  in  both  analyses.  Based 
on  the  analysis  results  shown  in  figure  D6,  it  demonstrates  that  the  simplified  crack 
qrowth  analysis  method  based  on  an  equivalent  stress  converted  from  the  stress 
spectrum  produces  the  same  crack  growth  life  as  that  using  the  cycle-by-cycle  method 

Suppose  the  actual  da/dN  data,  instead  of  that  from  the  Walker’s  equation,  were  used 
in  the  cycle-by-cycle  method.  The  crack  growth  life  obtained  is  compared  with  that  from 
the  simplified  method  in  Figure  7.  It  shows  that  the  crack  growth  life  obtained  from  the 
cycle-by-cycle  method  (144,757  flights)  is  about  3.5%  shorter  than  that  from  the 
simplified  method  (149,850  flights). 
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